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Magnetic order in the thermally quenched photomagnetic Prussian blue analogue coordination 
polymer Ko.27Co[Fe(CN) e ]o.73[D206]o.27-l-42D20 has been studied down to 4 K with unpolarized 
and polarized neutron powder diffraction as a function of applied magnetic field. Analysis of the 
data allows the onsite coherent magnetization of the Co and Fe spins to be established. Specifically, 
magnetic fields of 1 T and 4 T induce moments parallel to the applied field, and the sample behaves 
as a ferromagnet with a wandering axis. 

PACS numbers: 75.50.Xx, 75.25.-j, 75.30.Gw, 75.50.LK 



I. INTRODUCTION 

Manipulating magnetization with photons is now a ma- 
jor research focus because it may yield materials capable 
of dense information storage. An epitomic example of a 
photomagnetic coordination polymer is potassium cobalt 
hexacyanoferrate, K Q Co[Fe(CN)6]^-nH20 (from now on 
referred to as Co-Fe, with the crystal structure shown 
in Fig. [T]), which displays magnetic order and an op- 
tical charge transfer induced spin transition (CTIST). 1 
The details of the magnetism in Co-Fe have been inves- 
tigated with bulk probes such as magnetization,^ and 
AC-susceptibility, 3 as well as atomic level probes such 
as X-ray magnetic circular dichroism (XMCD)p' and 
muon spin relaxation (/z-SR)P However, we utilize neu- 
tron scattering because it is capable of extracting the 
magnetic structure, including the length and direction of 
the magnetic moments associated with different crystal- 
lographic positions. 

Neutron scattering research has been important in un- 
derstanding the structure of materials similar to Co- 
Fe. For example, neutron diffraction has been used 
to elucidate the location of water molecules, to iden- 
tify the long-range magnetic order, and to explore the 
spin delocalizetion in Prussian blueP^ Later work used 
similar techniques to investigate hydrogen adsorption 
in Cu3[Co(CN) 6 ] 2 , along with vibrational spectroscopy, 9 
and neutron vibrational spectroscopy was also measured 
in Zn3[Fe(CN) 6 ]2.^ Likewise, magnetic structure de- 
termination with neutrons was used to explore neg- 
ative magnetization in Cu . 73 Mn . 77 [Fe ( CN) 6 ] ■ nH 2 O 
and to extract on-site moments in Berlin greerP^ESl and 
in (Ni K Mni_ 2; )3[Cr(CN)6]2 molecule-based magnets.^ 

To this end, we have performed neutron powder diffrac- 
tion (NPD) on deuterated Co-Fe samples in magnetic 
states resulting from thermal quenching. Briefly, at room 



temperature, the photomagnetic Co-Fe with optimal iron 
vacancies is paramagnetic, with a transition to the dia- 
magnetic low-spin state when cooling below nominally 
200 K. It is below 100 K that applied light may convert 
molecules from diamagnetic to paramagnetic and back, 
and below around 20 K where this effect is most strik- 
ing due to the large susceptibility of the magnetically 
ordered state. However, a magnetically ordered state 
may also be ac hieved at low temperatures by thermally 
quenching, 1 5 1 1 where the paramagnetic 300 K state is 
cooled so quickly to below 100 K that it does not re- 




FIG. 1. (Color online) Co-Fe unit cell. Crystallographic posi- 
tions of atoms within the unit cell are illustrated with cubes 
for K (cyan), Co (red), Fe (blue), C (white), N (black), and 
coordinated D2O, O positions (yellow), while the interstitial 
D2O density is displayed using contoured isosurfaces (green) 
and Fe-C bonds are displayed as tubes (white). Details of 
structure determination are presented in Section III. 
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lax to the diamagnetic ground state. It is this magnetic, 
thermally quenched state that we study with NPD as a 
function of magnetic field, while complementary magneti- 
zation, transmission electron microscopy (TEM), Fourier 
transform infrared spectroscopy (FT-IR), and elemental 
analysis have also been performed on the sample. We 
find that Co-Fe possesses a correlated spin glass ground 
state that is driven via magnetic field to behave as a fer- 
romagnet with a wandering axis. 

II. EXPERIMENTAL DETAILS 33 

A. Synthesis 

To begin preparation of K Q Co[Fe(CN) 6 ] J g-nD 2 pow- 
der, a 75 mL solution of 0.1 mol/L KNO3 in D2O 
was added to a 75 mL solution of 20xlCP 3 mol/L 
K 3 [Fe(CN) 6 ] in D 2 0, and stirred for ten minutes. While 
continuing to stir, a 300 mL solution of 5xl0~ 3 mol/L 
C0CI2 in D2O was added drop- wise over the course of two 
hours. Stirring of the final solution was allowed to con- 
tinue two additional hours subsequent to complete mix- 
ing. Next, the precipitate was collected by centrifugation 
at 2000 rpm (210 rad/s) for 10 min (600 s) and dried un- 
der vacuum. This procedure was repeated 14 times until 
4.37 g of powder was collected. Potassium ferricyanide, 
anhydrous cobalt chloride, and potassium nitrate were 
all purchased from Sigma- Aldrich. To remove water, 
the potassium nitrate was heated in an oven to 110 °C 
(383 K) for 4 h before use. All other reagents and chem- 
icals were used without further purification. Deuterium 
oxide was purchased from Cambridge Isotope Laborato- 
ries, Inc. 



ference device (SQUID) magnetometer. Infrared spectra 
were recorded on a Thermo Scientific Nicolet 6700 spec- 
trometer. Energy dispersive X-ray spectroscopy (EDS) 
and TEM were conducted on a JEOL 2010F super probe 
by the Major Analytical Instrumentation Center at the 
University of Florida (UF). The UF Spectroscopic Ser- 
vices Laboratory performed combustion analysis. 



C. Analysis Preparations 

For NPD, 4.37 g of powder were mounted in a cylin- 
drical aluminum can. Thermal quenching to trap the 
magnetic state was achieved by filling the cryostat bath 
with liquid helium and directly inserting a sample stick 
from ambient temperature. To avoid hydrogen impuri- 
ties, the powder was wetted with deuterium oxide, and to 
avoid sample movement in magnetic fields, an aluminum 
plug was inserted above the sample. To measure mag- 
netization, samples heavier than 10 mg were mounted in 
gelcaps and held in plastic straws. Thermal quenching 
in the SQUID was achieved by equilibrating the cryo- 
stat to 100 K, and directly inserting the sample stick 
from ambient temperature. For measurements in 10 mT, 
samples are cooled through the ordering temperature in 
10 mT, and for 1 T and 4 T measurements, there is no 
observed thermal hysteresis. For FT-IR, less than 1 mg 
amounts of sample were suspended in an acetone solution 
and deposited on KBr salt plates and allowed to dry. For 
EDS and TEM, acetone suspensions of the powder were 
deposited onto 400 mesh copper grids with an ultrathin 
carbon film on a holey carbon support obtained from Ted 
Pella, Inc. 



B. Instrumentation 



D. Diffraction analysis scheme 



Neutron powder-diffraction experiments were con- 
ducted using the HB2A diffractomcter at the High Flux 
Isotope Reactor ,121 using the Ge[113] monochromator 
with A = 2.41 A (0.241 nm). Sample environment on 
HB2A utilized an Oxford 5 T, vertical-field magnet with 
helium cryogenics. Neutron polarization was achieved 
with a 3 He cell that produced 79% polarization at the 
beginning of the experiment and decayed to 63% polar- 
ization after 20 hours at the end of the experiment, to 
give an average polarization of 71% for both up and down 
polarization measurements, and we did not perform po- 
larization analysis after the sample but instead followed 
established m ethod s for powder diffraction with polar- 
ized neutrons! 19 * 20 * The flipping difference spectra were 
obtained by subtracting the diffraction data, measured 
with the incident neutron polarization parallel to the ap- 
plied field and magnetization, from the data recorded 
with the incident polarization antiparallel to the field. 
Magnetic measurements were performed using a Quan- 
tum Design MPMS XL superconducting quantum inter- 



Intensities were fit to the standard powder diffraction 
equation with a correction for absorption, 



v ' sin6» sin26» n ; 



with 

V (0) = e~ 



(1) 



,(2) 



where A is an overall scale factor, m^ki is the multiplic- 
ity of the scattering vector, F is the structure factor, 6 
is the scattering angle, fi is the linear attenuation coef- 
ficient, and R is the radius of the sample cylinder. For 
our sample and experimental arrangement, fiR = 0.17, 
which has little effect on the observed intensities aside 
from scale. The structure factor has nuclear (-Fjv) and 
magnetic (Fm) contributions, and for unpolarized neu- 
trons 
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On the other hand, Fn and Fm can coherently interfere 
for polarized neutrons such that, for moments co-linear 
with P, 

\F\ 2 = \F N +F M \ 2 = \F N \ 2 + \F M \ 2 ±2PF N F M , (4) 

where P is the neutron polarization fraction and the 
sign of the final term depends upon up or down neutron 
polarizationPH For nuclear scattering, 

F N (hkl) = Y,".; 1 ',''"' 1 ' U . (5) 

3 

where the sum is over all atoms in the unit cell, n is 
related to the average occupancy, b is the coherent nu- 
clear scattering length, G is the hkl dependent recipro- 
cal lattice vector, d is the direct space atomic position, 
and W = BQ 2 /16tt 2 is the Debye- Waller factor. For 
magnetic scattering, all coherent scattering is modeled 
to be along the applied field, which is perpendicular to 
the scattering plane, so that 

F M (hM) = i^Y^ m i(Qy G ' di t~ Wi > w 

3 

where = 2.695 fm, and the magnetization can be 
written as 

m,(Q) = (LjjfaiQ) + 2{S z ) 3 fsAQ) 
= 9J,j( J z)jfj,j(Q) 

= (JMaLjhAQ) + gsjfsAQ)) > (7) 

where (J z ) is the average total angular momentum, (L z ) 
is the average orbital angular momentum, (S z ) is the av- 
erage spin angular momentum, fj(Q) is the magnetic 
form factor for the total angular momentum, is 
the magnetic form factor for the orbital angular momen- 
tum, fs(Q) is the magnetic form factor for the spin an- 
gular momentum, and oj, g^ and gs may be determined 
by Wigner's formulaP^l The tabulated form factor values 
within the dipole approximation are used for the spin 
and orbital form factorsPS Squared differences between 
observed and calculated intensities were minimized using 
a Nelder-Mead simplex algorithm. Open-source Python 
2.7 libraries were utilized to aid in plotting routines, 
matplotlib 1.0.1 and Mayavi2, and computation, NumPy 
1.6.1 and SciPy 0.7.2. Reported uncertainties of fit pa- 
rameters are the square root of the diagonal terms in the 
covariance matrix multiplied by the standard deviation 
of the residuals. 



III. RESULTS AND ANALYSES 

The nuclear crystal structure of Co-Fe can be mod- 
eled with space group Fm3m (No. 225), where ferri- 
cyanide molecules and cobalt ions are alternately cen- 
tered on the high symmetry points of the unit cell, with 
heavy- water bound to cobalt when ferricyanide is absent, 
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FIG. 2. (Color online) Neutron powder diffraction of Co-Fe 
at T = 40 K. Observed scattering is shown as open circles 
(obs), a full fit including sample mount contributions is shown 
as a black line (calc'd), the diffuse background is illustrated 
with a magenta line (bgr), and the residuals of the fit are 
shown below the zero line with a green line (resid). The signal 
due to Co-Fe is emphasized with a red filling. Experimental 
uncertainties derived from counting statistics are smaller than 
the plotting symbols. 



TABLE I. Atomic coordinates and occupancies for Co-Fe at 
T = 40 K. 



atom 


position 


n 


X 


y 


z 


Co 


4a 


1 


0.5 


0.5 


0.5 


Fe 


4b 


0.73 











C 


24e 


0.73 


0.212 








N 


24e 


0.73 


0.313 








K 


8c 


0.135 


0.25 


0.25 


0.25 


O 


24e 


0.27 


0.243 








D 


96k 


0.135 


0.303 


0.060 


0.060 



and potass ium ions and heavy-water molecules filling in 
voids! 6 ! 8 * 24 ! This structure is used as a starting point to 
fit the T — 40 K thermally quenched Co-Fe contribution 
to the measured intensity profile, Fig. [2j which also has 
sample mount contributions due to P63/ mm c (No. 194) 
D 2 and Fmim (No. 225) aluminum! 25 1 26 1 Incomplete 
trapping of the high-temperature state in Co-Fe gives 
rise to a highly microstrained nuclear structure,^ and 
to account for this effect during refinement, we use an 
asymmetric double sigmoidal peak shape, namely 

I f _ 1 \ / 1 \ 
Va2s ~ 2T9^ [} i + e-^t) [l + e-^J ' 

where I is the intensity, w is the width, and 9 C is the 
center of the reflection, and these fits yield an effective 
lattice constant of 10.23 A. Observed Co-Fe reflections 
that can be clearly separated from sample holder reflec- 
tions are used to extract structure factors. In modeling 
the unit cell, the cobalt to iron ratio was determined 
with EDS, while the room temperature oxidation states 
with FT-IR. The carbon and nitrogen content were es- 
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FIG. 3. (Color online) Magnetic neutron powder diffraction of Co-Fe at T — 4 K as a function of applied magnetic field. The 
difference between the T — 40 K diffractogram and the T = 4 K diffractogram (open circles) along with profile fits to intensities 
(red line from model #1 for 4 T and 1 T data in Table II) are shown for (a) 4 T, (b) I T, and (c) 10 mT. Additionally, the 
difference between the T = 4 K up-neutron-polarization diffractogram and the T — 4K down-neutron-polarization diffractogram 
(open circles) along with profile fits to intensity (red line from model #5 for polarized data in Table II) are shown for (d) 1 T. 
Uncertainty bars on experimental data points are statistical in nature representing one standard deviation from the mean, 
using counting statistics. 



tablished with combustion analysis, and the potassium 
ions provide charge balance. Finally, the heavy- water 
concentration and positions were refined along with the 
scale factor to fit the structure factors. 

To begin, refinement yielded interstitial heavy-water 
pseudo-atoms at the 8c position (n = 0.618, B — 
5) and the 32f position (x = 0.3064, n = 0.333, B 
= 20), after which all other parameters were fixed 
(Table [I]) and Fourier components of the heavy-water 
were further refined to give the interstitial distribu- 
tion shown in Fig. [I] These refinements give a chemi- 
cal formula of Ko.27Co[Fe(CN) 6 ]o.73[D 2 6 ]o.27-l-42D 2 0. 
Moreover, at room temperature, the more com- 
plete chemical formula with oxidation states of the 
metal ions included is K .27Coo1,4CoQ^ 6 [Fe 3+ (CN) 6 ]o.58 
[Fe 2+ (CN) 6 ] .i5[D 2 O 6 ] .27-l-42D2O, or more compactly 
represented by Coo^Co^gFe^gFe 2 /^. 

Having highly ionic wavefunctions, the magnetic 
ground states of Co and Fe in Co-Fe are well de- 
scribed with ligand field theoryP^ As displayed in Fig. [l] 
the iron atoms are octahedrally coordinated by carbon 
atoms that introduce a ligand field splitting parame- 
ter (Ap e ) of approximately 0.70 aJ (35,000 cm -1 or 4.3 
eV), and typical Fe Racah parameters put d 5 — Fe 3+ into 
a 2 T2 g ground state, and d 6 — Fe 2+ into a diamagnetic 



A\ g ground state. Similarly, cobalt atoms arc octa- 
hedrally coordinated with oxygen and nitrogen atoms 
to give A Co 3+ « 0.46 aJ (23,000 cm" 1 or 2.9 eV) 
for d 6 — Co 3+ that has a diamagnetic 1 A\ g ground state, 
and A Co 2+ w 0.20 aJ (10,000 cm" 1 or 1.2 eV) for 
d 7 — Co 2+ that has a 4 T iff ground state, using typical 
Co Racah parameters. At temperatures much less than 
the spin-orbit coupling energy, only the lowest energy 
total angular momentum levels are appreciably occu- 
pied, so that the relevant states are Fe 3+ [J = 1/2, 
.gj = (2+4fc)/3, g L = 4fc/3, gs = 2/3] and Co 2 +[J = 1/2, 
.9,7 = (10 + 2Afe)/3, g L = 2Afc/3, g s = 10/3], where 
A is expected to be nearly 1.5 due to the weak ligand 
field, and k is the orbital reduction parameter. It is 
worth noting that analogous orbitally degenerate terms 
have been observed for d 5 -Fe 3+ ( 2 T 2g ) in K 3 Fe(CN) 6 pSl 
and d 7 -Co 2 +( 4 T lg ) in K 1 . 88 Co[Fe(CN) 6 ] .g 7 -3.8H 2 O and 
Nai. 5 2K .o4Co[Fe(CN) 6 ]o.89-3.9H20P Alternatively, if 
interaction with the lattice drastically quenches the or- 
bital moment, spin-orbit coupling no longer splits the 
ground states and the magnetic parameters become 
Fe 3 +[J = 1/2, 5 ,7 = 2, g L = 0, g s = 2] and Co 2 +[J = 
3/2, gj = 2, gL = 0, gs = 2]. To estimate the rel- 
ative proportion of the different oxidation states in the 
thermally quenched state, the effective paramagnetic mo- 
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TABLE II. Comparison of the eight magnetic models, as described in the text, numbered (#) 1 — 8 for Co-Fe at T — 4 K 
in different magnetic fields tabulated as "cond.", which is shorthand for experimental condition, where "P" designates the 
data acquired with polarized neutrons. Here, "align." is short for "moment alignment," where + denotes parallel alignment 
of moments and - denotes antiparallel alignment of moments. The units of m Zt co and m Zi p e are /is, and the units of M are 
\ib mol -1 . The sum of the residuals are normalized to model 1 for each experimental condition. 



# 


cond. 


align. 


9s, Co 


9l,Co 






J z Co 




77l z q 


Tfl z F & 


M 


residual 2 


1 


4 T 


4- 


10/3 


1 


2/3 


4/3 




1 3 + 03 


2.7 


0.3 


2.6 


1.000 


2 


4 T 




10/3 


1 


2/3 


4/3 


0.62 ± 0.02 


0.00 ± 0.02 


2.7 


0.0 


2.4 


1.019 


3 


4 T 


_|_ 


10/3 


1 


2 


o 


64 + 02 


1 2 + 04 

V 7 . 1 _ 1 VI . VII 


2.8 


0.2 


2.6 


1.000 


4 


4 T 




10/3 


1 


2 





0.63 ± 0.02 


0.00 ± 0.11 


2.7 


0.0 


2.5 


1.020 


5 


4 T 


4- 


2 


o 


2/3 


4/3 


1 08 + 03 


16 + 03 


2.2 


0.3 


2.1 


1.028 


6 


4 T 




2 





2/3 


4/3 


1.04 ± 0.03 


0.00 ± 0.03 


2.1 


0.0 


1.9 


1.049 


7 


4 T 


_|_ 


2 


o 


2 


o 


1 OQ + 03 


1 2 + 04 


2.2 


0.2 


2.1 


1.030 


g 


4 T 




2 


o 


2 


n 


i nc _i_ n 03 


o no + n 03 


2.1 


0.0 


1.9 


1.048 


1 


1 T 


4- 


10/3 


1 


2/3 


4/3 


37 + 03 


20 + 0Q 


1.6 


0.4 


1.7 


1.000 


2 


1 T 




10/3 


1 


2/3 


4/3 


0.37 ± 0.02 


0.00 ± 0.04 


1.6 


0.0 


1.4 


1.020 


3 


1 T 


_|_ 


10/3 


1 


2 


o 


38 + 03 


1 2 + 07 


1.6 


0.2 


1.6 


1.004 


4 


1 T 




10/3 


1 


2 





0.38 ± 0.05 


0.00 ± 0.11 


1.6 


0.0 


1.5 


1.019 


5 


1 T 


i 


2 


o 


2/3 


4/3 


64 + 04 


20 + OQ 


1.3 


0.4 


1.4 


1.001 


6 


1 T 




2 





2/3 


4/3 


0.64 ± 0.10 


0.00 ± 0.21 


1.3 


0.0 


1.2 


1.022 


7 


1 T 


_i_ 


2 


o 


2 


o 




1 3 + 07 

yj . id _l_ U . \J i 


1.3 


0.3 


1.3 


1.006 


8 


1 T 




2 





2 





A GO _l_ A A A 

0.63 ± 0.09 


A AA _l_ A AO 

0.00 ± 0.08 


1.3 


0.0 


1.1 


1.021 


1 


P 


+ 


10/3 


1 


2/3 


4/3 


0.38 ± 0.01 


0.06 ± 0.02 


1.6 


0.1 


1.5 


1.000 


2 


P 




10/3 


1 


2/3 


4/3 


0.40 ± 0.02 


0.00 ± 0.01 


1.7 


0.0 


1.5 


1.058 


3 


P 


+ 


10/3 


1 


2 





0.36 ± 0.02 


0.18 ± 0.05 


1.6 


0.4 


1.6 


1.004 


4 


P 




10/3 


1 


2 





0.40 ± 0.02 


0.00 ± 0.02 


1.7 


0.0 


1.5 


1.065 


5 


P 


+ 


2 





2/3 


4/3 


0.68 ± 0.02 


0.06 ± 0.02 


1.4 


0.1 


1.3 


0.995 


6 


P 




2 





2/3 


4/3 


0.71 ± 0.02 


0.00 ± 0.02 


1.4 


0.0 


1.3 


1.058 


7 


P 


+ 


2 





2 





0.66 ± 0.03 


0.16 ± 0.06 


1.3 


0.3 


1.3 


1.005 


8 


P 




2 





2 





0.74 ± 0.03 


0.00 ± 0.05 


1.5 


0.0 


1.3 


1.079 



merit is linearized as a function of temperature for the 
300 K and 100 K states,^ and the measured lattice con- 
stant is compared to a weighted average of quenched and 
ground-state lattice constants,^ to self-consistently give 
( -' ol30^ o/io^ e o^54^ e o/i9 f° r tne magnetic quenched state 
at 100 K and below that is analyzed in detail herein. 

Cooling the sample further, subsequent to quenching, 
the bulk magnetization measured in 10 mT showed the 
well-documented upturn at around 15 K corresponding 
to the onset of magnetic order. Therefore, additional 
NPD was performed at 4 K in applied fields of 10 mT, 
1 T, and 4 T, Fig. |3j to compare to the scattering in the 
paramagnetic state. Furthermore, polarized NPD was 
performed at 4 K in an applied held of 1 T, Fig. [3] (d), 
where the difference between diffractograms for up and 
down incident neutron polarizations increases signal to 
noise of the measured magnetic structure at reflections 
with large nuclear contributions. For each of the three 
experimental conditions where magnetic scattering is ob- 



served, we compare the results of eight plausible but dif- 
ferent models that all have moments along the applied 
held. Specifically, each possible case considers various 
combinations of the parallel or antiparallel alignment of 
Co and Fe moments when each ion possesses either spin- 
only or orbitally degenerate magnetic states, Table |TT] 
The analyses indicate that most magnetism resides on 
the Co 4a site for all models, with a parallel alignment of 
Fe and Co moments giving the best hts and x 2 surfaces 
suggesting a reduced but present orbital moment on both 
ions. No magnetic scattering is observed in 10 mT, and 
increased coherent magnetic scattering appears with in- 
creasing field, which is consistent with the presence of sig- 
nificant random anisotropy, where a correlated spin glass 
(CSG) is the ground state and sufficiently large fields 
cause entrance into a ferromagnetic phase with wander- 
ing axis (FWA) state or at even larger helds a nearly 
collinear (NC) stateP^ 

Analytical expressions for the magnetization process 
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for magnets with random anisotropy are available^ for 
the Hamiltonian 

-D C J2(SD 2 - 9HbY, H ■ & > ( 9 ) 

i i 

where J is the superexchange constant, S is the spin 
operator, D r is strength of the random anisotropy, h is 
the direction of the random on-site anisotropy, D c is the 
strength of the coherent anisotropy, g is the Lande factor, 
and H is the applied field. In the FWA regime where the 
applied field energy is larger than the random anisotropy 
energy but much less than the exchange field, the low 
temperature magnetization is 



Mfwa = M s 



Ms 



QV2D 2 ttM s 



5<K 2 a 3 (zJ) 3 / 2 (H + H c y/ 2 



1 - 



D 



1/2 
FWA 



{H + He) 1 / 2 



(10) 



where Ms is the saturation magnetization, f2 is the 
integrated local anisotropy correlation function, z is 
the number of magnetic neighbors, a is the mean dis- 
tance between neighboring spin sites, He is the coher- 
ent anisotropy field, and Dp wa is a measure of the ran- 
dom anisotropy to superexchange strengths. For the NC 
phase that is reached when the applied field and coherent 
anisotropy field are much larger than the exchange field, 
the low temperature magnetization is 



M 



NC 



Ms 



M, 



AD^Ms 



15a 6 (H 
1 

(H + He) 2 



He) 2 



D 2 

^NC 



(11) 



where D^e is a measure of the random anisotropy. Based 
upon the magnetic ordering temperature, Co-Fe is ex- 
pected to be in an FWA-like phase, and although both 
FWA and NC expressions may be fit to the low tempera- 
ture magnetization data, Fig.[4j the parameters extracted 
by the NC fit are not consistent with the derivation limit 
for Eq. 11, further suggesting an FWA-like state. 




SQUID magnetization 
■ model fit 



0.0 



FWA: ^ = 2.7,0^=1.0,/^ =1.6 
NC: M = 2.1, D.„ = 5.2, H„ = 6.4 
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FIG. 4. (Color online) SQUID magnetization of Co-Fe at 
T — 2 K. Experimental magnetization is shown as open cir- 
cles (SQUID magnetization) and model fits as a red line 
(model fit), where FWA and NC are visually indistinguish- 
able; M(l T) = 1.1 ± 0.1 fi B moU 1 and M(4 T) = 1.6 ± 
0.2 ixb moP 1 . The units of Ms are hb mol -1 , and He, 
Dfwa, and Dnc are all units of Tesla. Uncertainty bars rep- 
resent one standard deviation from the mean, where statistics 
are generated by measuring the magnetization of the 14 syn- 
thesis batches required to generate 4.37 g for NPD. 
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IV. DISCUSSION 

We have presented neutron diffraction 
and bulk magnetization measurements of 
Ko.27Co[Fe(CN) 6 ]o. 73 [D 2 6 ]o.27-1.42D 2 that sug- 
gest a CSG ground state that enters a FWA-like state 
in applied magnetic fields of the order 1 T and larger 
(Fig. [5]). This conclusion is based upon the field 
dependence of the magnetization, and particularly the 
diffraction experiments that show an absence of long 
range order in the 10 mT data and ordered moments 
that are induced by the applied field, ruling out a 



FIG. 5. (Color online) An illustration of magnetic structure 
for different magnetic field regimes. Here, the magnetic field 
points towards the top of the page, short arrows represent iron 
moments, and long arrows represent cobalt moments, (a) The 
Co-Fe sample cooled in zero field has a CSG-like state with 
no average on-site moment, as shown for the measurement of 
magnetic scattering in 10 mT, Fig. 3 (c). (b) The application 
of magnetic field cants the moments towards the field (Fig. 3 
(b) and (d) ), and (c) larger fields induce larger average mo- 
ments (Fig. 3 (a) ). (d) More complicated mesoscopic states 
that contain texture are also possible, but are not unambigu- 
ously determined with our data. 
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high-field domain magnetization process. This random- 
anisotropy-based magnetization process explains the 
appreciable slope observed for Co-Fe even in fields of 7 
T at 2 K (Fig. 4) , in a way similar to the magnetization 
process at low temperature and high magnetic fields 
for the vanadium tetracyanoethylene molecular magnet 
prepared using solvent based methods.^ This mag- 
netization process is different than for other reported 
cubic complex cyanide systems that have magnetically 
ordered ground states and saturate magnetization at 
2 K and 7 TpHUl The CSG ground state is consis- 
tent with previous AC-susceptibility measurements of 
Ki_2xCoi +!B [Be(CN)a]-vHaO (0.2 < x < 0.4, y « 5) 
that showed glassy behavior,^ although the relative 
orientation we find for Co and Fe at 1 T is contrary 
to the XMCD experiment that reported antiparallel 
Co and Fe at 1 T in Rbi. 8 Co 4 [Fe(CN) 6 ] 3 . 3 -13H 2 and 
K .iCo 4 [Fe(CN) 6 ] 2 . 7 -18H 2 o! 31 For fields of 1 T and 4 T, 
we find a parallel alignment of Co and Fe moments 
minimizes the residuals between model and data, and 
this alignment is clearly seen for the low-angle 4 T peaks, 
Fig. [6] However, the lack of coherent scattering for 
CSG means we cannot strictly discuss the nature of the 
superexchange interaction in the ground state, although 
we infer some significant ferromagnetic character based 
upon the high-field regime. One must also be careful 
about applying qualitative Goodenough-Kanamori rules 
to this systempSlEU as the single-electron states are not 
only mixed from electrostatic interactions, but also due 
to the aforementioned presence of spin-orbit coupling. 
A band structure calculation using a full potential 
linearized augmented plane wave method resulted in 
an antiferromagnetic ground state for Co-Fe, with 
+0.296 [i B on the Co site and -0.280 /ig on the Fe sitej^ 
although such values do not agree with experimental 
findings. 

For the Co-Fe presented, care was taken to ensure 
that the average particle size was greater than sslOO nm 
(as measured by TEM) to avoid finite size effects,^ 
and the FWA-like phase can explain the previously 
reported changes in low temperature high field mag- 
netization with particle size as a tuning of the local 
random anisotropy with size, larger particles requiring 
higher fields to saturate at base-temperature as they are 
deeper in a glassy phase. The saturation value for bulk 
K . 27 Co[Fe(CN) 6 ] . 73 [D 2 6 ]o.27 • 1.42D a O, M s = 2.7 ± 
0.3 \ib mol -1 , is comparable to a variety of similar states 
with different degrees of orbital reduction on Co and Fe 
sites; exempli gratia, considering complete orbital mo- 
ments on both Co and Fe gives Ms ~ 2.5 /is mol -1 and 
spin-only moments gives Ms ~ 3.2 /ig mol -1 . 

The NPD experiments show a model-independent in- 
crease in coherent magnetic scattering as a function of 
field, and a slightly form-factor dependent ratio of Co to 
Fe moments, Table |H| A parallel alignment of Co and Fe 
is found for both 1 T and 4 T, and when an antiparallel 
alignment is forced, Fe moments go to zero to achieve best 
fits. The moment ratio is heavily dictated by the low-Q 




22 23 24 25 26 27 28 29 
20 (degrees) 



FIG. 6. (Color online) Visual comparison of magnetic con- 
figurations for Co-Fe. Here, the low angle H = 4 T data 
(open circles) are set side by side with model #1 for 4 T from 
Table II (thick-solid, red line), model #2 for 4 T from Ta- 
ble II (dotted, black line), a ferrimagnetic structure with a 
3:1 Co:Fe spin ratio (dashed, blue line), and a ferromagnetic 
structure with a 3:1 Co:Fe spin ratio (thin-solid, black line). 
Uncertainty bars are representative of one standard deviation 
from the mean, using counting statistics. 



peaks where the form- factor has little effect, while the 
scale of the moments is different depending upon the pre- 
sumed shape of the scatterer. Previous neutron diffrac- 
tion measurements have shown covalency effects, due to 
CT-bonding and 7r-back-bonding with CN, to be important 
in the chemically similar CsK 2 [Fe(CN)6].^ Covalency 
can increase the direct-space size of the moments, thereby 
decreasing the reciprocal-space size even in the presence 
of orbital moments. For Co-Fe, the smaller reciprocal- 
space form-factors give magnetizations most similar to 
those determined by SQUID, although covalency makes 
assignment of orbital and spin magnetism based upon 
spatial distribution inconclusive. We do not refine the 
form-factor in this manuscript because high parameter 
covariance is introduced. Finally, small quantitative dif- 
ferences between SQUID and NPD moment values may 
also be due to sample inhomogeneity overestimating the 
NPD moments, and unitemized experimental uncertain- 
ties due to the complicated and highly un-stoichiometric 
formulation of the Co-Fe material, but our conclusions 
remain robust with respect to such perturbations. The 
line-widths for magnetic and nuclear NPD are similar, 
suggesting comparable domain sizes for the scattering ob- 
jects, but it is possible that the induced moments have a 
texture over some other length scale, Fig. 5 (d), a possi- 
bility suggested by cluster-glass behavior in AC suscep- 
tibility studiesP As shown in the Appendix, regions of 
coherent magnetization at an angle £ from the applied 
field would rescale the measured NPD longitudinal mo- 
ment by c° s £^/~: os 2 2 g +1 , but the similarity between the 
polarized and unpolarized magnetic diffraction further 
suggests that such an effect is small in our samples. 



V. CONCLUSIONS 

The neutron diffraction and bulk magnetization mea- 
surements of Co-Fe suggest a magnetization process 
that evolves from a correlated spin glass to a quasi- 
ferromagnetic state with increasing magnetic field, where 
average Co and Fe moments are induced to lie along the 
applied field. When considering memory storage appli- 
cations of molecule based magnetic materials, structure- 
property relationships that may give rise to coherent and 
random anisotropy will be important to consider. 
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Appendix: Effect of Canting on Intensity 

A powder sample consisting of domains canted at an 
angle £ away from the applied field, with random rota- 
tional distribution, may give the same unpolarized NPD 



signal as domains along the field, but with a different 
magnetic moment. With the magnetic field along the 
z-axis, and the scattering vector along the x-axis, the 
uncanted magnetic moment is simply 

M u = (0,0, M) , (A.l) 

and the canted moment can be expressed as 

M can ted = M(sin £ cos </>, sin £ sin </>, cos £) , (A. 2) 

where £ is the canting angle, <f> is the rotation angle about 
the field, and M is the magnitude of the magnetic mo- 
ment. The interaction vector^ 

l M ^! 2 = ^(8ap-Q a Qp)MZMi) , (A.3) 

a,f) 

then gives a dependence of the intensity on the canting 
angle such that 

|M U , ± | 2 = M 2 , (A.4) 

and 

|M canted ,J 2 = M 2 ^!|±l , (A .5) 

where random 0-angles have been averaged over. There- 
fore, an uncanted model with a z-component (relevant to 
compare with longitudinal magnetization) of M can give 
the same unpolarized NPD intensity as a canted model 

with a z-component of Mcos£* / cos 2 2 f+1 - 
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